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Abstract. We analyze the dynamics of up to 10 5 electrons resulting from illuminating a Xenon cluster with 
9093 atoms with intense laser pulses of different length and peak intensity. Interesting details of electron 
motion are identified which can be probed with a time resolution of 100 attoseconds. Corresponding 
experiments would shed light on unexplored territory in complex electronic systems such as clusters and 
they would also permit to critically access the present theoretical description of this dynamics. 

PACS. 36.40.Gk Plasma and collective effects in clusters - 31.15.Qg Molecular dynamics and other nu- 
merical methods - 36.40.Wa Charged clusters - 33.80.Wz Other multiphoton processes 
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1 Introduction 

Recent advance in laser technology has led to the creation 
of subfemtosecond (100 attosecond) strong laser pulses 
which allow one in the future to resolve much shorter time 
scales in microscopic dynamics than it has been pos- 
sible so far. More specifically, while until know vibrational 
motion in molecules could be resolved in time experimen- 
tally, it will be possible to follow electronic motion with 
attosecond resolution. Still, the atomic time unit, i.e., the 
period of the ground state electron in hydrogen, is of the 
order of 10 -17 s, i.e., roughly 10 attoseconds while for now 
a pulse length of 100 attoseconds and more is experimen- 
tally feasible. However, there are alternative methods such 
as tomographic imaging which use the high harmonic spec- 
trum of a femtosecond pulse to get attosecond resolved 
dynamics. In this way it could be demonstrated that the 
electronic wavefunction of a molecule could be directly im- 
aged This is fantastic progress, yet one might argue 
that we know in principle the electronic wavefunction of 
a diatomic molecule. We can even calculate it, and as it 
turns out experimental imaging and theoretical calcula- 
tions via a standard solution of the Schrodinger equation 
agree well [3]. 

On the other hand, for systems more complex than 
a diatomic molecule, e.g., a cluster consisting of many 
atoms, we neither know with certainty theoretically the 
dynamics of the electrons (we have to make model as- 
sumptions) nor do we have a way so far to access the 
electron dynamics time resolved in the experiment. 

Pioneering experiments 4 have demonstrated non-tri- 
vial time dependent dynamics in cluster motion under a 
strong laser pulse. This interesting dynamics can be at- 
tributed to resonance like behavior where either the next 
neighbor ions contribute in a cooperative way to enhance 



ionization |2j or the electrons, still bound to the cluster but 
not to individual atoms, enter a collective phase of motion 
which is susceptible to resonant absorption of radiation 
[Sj. Besides that, a nonlinear excitation of this resonance 
could be relevant [7|. A slightly different situation results 
from irradiation with VUV pulses as realized with the free 
electron laser [jy. There, inverse Bremsstrahlung seems 
to be the main mechanism of energy absorption. Details, 
however, remain controversial [QlllOj . Detailed knowledge 
would require the experimental ability to time resolve the 
electron motion. Since the electrons are no longer strongly 
bound to a single ion but to the extended cluster, there 
typical time scale of motion is a bit slower which would 
be ideal for probing it with subfemtosecond laser light. 

It is the purpose of this paper to explore which kind of 
details the electron dynamics exhibits during a standard 
femtosecond laser pulse a length of the order of 100 fs. 
Specifically, we will concentrate on energy spectra of elec- 
trons from a Xenon cluster with 9093 atoms during the 
irradiation with a laser pulse of length from 25 to 400 fs 
and peak intensity ranging from 10 to 10 16 W/cm 2 . 



2 Theoretical method 

Before we start to discuss the spectra we give a brief ac- 
count of the theoretical method how these spectra have 
been obtained. The initial cluster configuration is derived 
assuming an icosahedral symmetry of the atoms. The faces 
of the nth shell contain (n+l)(n+2)/2 atoms. Keeping 
only those atoms which are inside a sphere of radius Rq = 
50 A yields a cluster with 9093 atoms. This configuration 
is relaxed using pairwise Lennard-Jones potentials to find 
the optimum interatomic distances thereby freezing the 
cluster geometry. 
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The laser pulse has the form of a Gaussian (we use 
atomic units if not stated otherwise) 



F(t) = z F t cos cot 



= zVl 



exp 



log 2 (a/I? 



cos cot 



(1) 



with half width T which measures the pulse length, peak 
intensity / and the laser frequency of lo — 0.058 a.u. (cor- 
responding to 780 nm wavelength). The light is linearly 
polarized along the z-axis, denoted by the unit vector z. 
To make the simulation tractable we only model explicitly 
the outermost bound electron of an atom/ion [Hj. In this 
way more and more electrons "are generated" due to the 
removal from their mother ion during the laser pulse. The 
mutual attractive Coulomb potentials have soft cores. All 
existing charged particles, ions and electrons, are propa- 
gated classically exposed to their mutual attraction and 
repulsion as well as to the coupling with the light pulse. 
Handling up to about ten electrons per atom implies a 
total of up to 10 5 particles to be propagated over rel- 
ative long times (of the order of hundred optical cycles) . 
This is not possible with standard molecular dynamics but 
requires tree-code techniques Alternatively, one can 
start from the particle-in-cell concept to handle clusters 
of this size \12\ . 



3 Energy spectra of cluster electrons during 
the laser pulse 

In this section we present and interpret the energy spec- 
tra of the electrons. An overview of the spectra is given 
in Fig.Q the parameters for each panel of figure Fig.Q]are 
provided in Tabled Basically, for each pulse length T the 
peak intensity increases by a factor of 5 from top to bot- 
tom (set A to C). For increasing pulse lengths T from left 
to right the fluency of the pulse remains the same. Hence, 
from one figure to the next the peak intensity decreases 
by a factor of 2 since the pulse length increases by a fac- 
tor of 2. All figures show distributions of total energies 
of the electrons on the y-axis ranging from — lOkeV to 
+5 keV with the total energy of electron i at position n 
with momentum pi defined as 
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with r a and q a the position and the charge of ion a, re- 
spectively. 

The time evolution of the spectrum is shown from — T 
to max{+T, 100 fs}. Of course, the calculation of the dy- 
namics sets in long before —T, when the pulse intensity 
is still negligible. It should be pointed out that Fig.Q] con- 
tains for each plot and at each time the energies of up 
to 10 5 electrons. Just the generation of such a plot in a 
reasonable time requires to use the tree-code information. 



Table 1. Laser pulse parameters for the 15 plots in Fig.0 
The peak intensity / given in units of 10 15 W/cm 2 and the 
pulse length T in fs enters Eq. Q. Furthermore, the quiver 
amplitude x^ — ^/i/ui 2 is given in A and the ponderomotive 
energy U p = J/4oj 2 in units of eV, whereby ui — 0.058 a.u., the 
laser frequency corresponding to 780 nm wavelength. 
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25 


50 


100 


200 


400 


set 


I 


3.2 


1.6 


0.8 


0.4 


0.2 


A 




46.7 


33.1 


23.4 


16.5 


11.7 




u P 


181 


90.7 


45.4 


22.7 


11.3 


set 


I 


16 


8 


4 


2 


1 


B 




105 


73.9 


52.3 


37.0 


26.1 




u P 


907 


454 


227 


113 


56.7 


set 


I 


80 


40 


20 


10 


5 


C 


x w 


234 


165 


117 


82.6 


58.4 




u P 


4536 


2268 


1134 


567 


284 



3.1 Common features in all spectra 

A first orientation in the plots reveals that the lowest 
electron energies go always through a minimum after the 
laser pulse has reached its maximum intensity. The mini- 
mum assumes quite different values, for the parameters 
in Fig.^ roughly between —2 and — 15keV. The nega- 
tive electron energies are due to the positive background 
charge Q t = J2a la of the cluster which is a direct con- 
sequence of ionized electrons which have left the cluster. 
They appear at positive energies in the plots. As stronger 
the cluster is ionized as deeper the remaining electrons 
are bound by the back ground charge. Upon explosion 
of the ions the background charge spreads out and con- 
sequently, the energy of the trapped electrons decreases 
for longer times. Evident are two preferred energy regions 
(with high intensities) for the electrons: Ionized electrons 
have excess energies close to E — 0. The trapped electrons 
prefer energies close to the lowest possible values. 



3.2 General trends as a function of peak intensity: 
The role of the quiver amplitude 

The population of the two preferred regions is most promi- 
nent in the lower left corner of Fig.^ i.e., for maximum 
peak intensity. For minimum peak intensity (upper right 
corner) the opposite trend can be identified: A substantial 
number of electrons fills the area (of negative energy) be- 
tween the two preferred regions. If one compares the equi- 
librium size of the cluster given by the radius Rq = 50 A 
with the quiver amplitude of an electron in the laser field 

— Vl/u 2 , listed in Tabled the reason for this trend 
becomes obvious: For largest intensities (lower left corner) 

> Rq and the ionized electrons are driven far outside 
the cluster. The opposite is true for low intensities (up- 
per right corner) with the lowest intensities where despite 
quivering in the laser field the electrons remain well inside 
the cluster, leading to a distribution of negative electron 
energies due to the attractive background charge. 
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Fig. 1. Energy distribution of electrons according to Eq. under illumination of a Xegog3 cluster with laser pulses of different 
pulse length T and peak intensity I as detailed in Tabled The pulse shape is given in Eq. 0. The visible area has range of 
energy on the y-axis between — lOkeV and +5keV, and a range of time on the x-axis between —T and max{T, 100 fs}. 
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Fig. 2. Phase difference between laser cycle and periodic response of electrons inside the cluster volume for the same laser pulse 
parameters as shown in Fig.0 



3.3 The trace of escaping single ions 

In the panels corresponding to large pulse lengths and 
high peak intensity (lower right corner of Fig.^) one sees 
individual lines of electron energies above the energy mini- 
mum and converging to each other for longer times. These 
lines are energies of electron being trapped in excited or- 
bits around a single ion due to rising interatomic barriers 
when the cluster expands. Hence, large pulse lengths are 
necessary for this feature. It also occurs in the upper right 
corner of Fig. ^ (i.e., for lower intensities) but it is masked 
there by other electrons with similar energies. They can- 
not leave the cluster since the laser intensity is too low. 



3.4 Atomic versus cluster effects 

Very short pulses in small clusters explore primarily atomic 
properties of the atoms and ions within the cluster since 
there is not enough time for a cooperative response of the 
cluster as a whole, e.g, by expansion or thermalized col- 
lective electron motion. Large clusters, however, develop 
upon the first atomic field ionization for each atom a sub- 
stantial positive background charge (in our case roughly 



Q = +10 4 for single ionization of the atoms). The main 
mechanism for further ionization is therefore the escape 
from the potential generated by the background charge. 
It can be modelled by assuming a homogeneously charged 
sphere with total charge Qt and (cluster) radius Rt, which 
leads to the potential 



Vt(r) 



-Qt{m 2 -r 2 )/(2R t 3 ) tfr<Rt 



-Qt/r 



if r > R t , 



(3) 



where quantities with subscript t are weakly time depen- 
dent. All electrons in this potential can leave the cluster 
if the height of the potential barrier at maximum field 
strength within a laser period is lower or equal to the po- 
tential minimum. This happens at the field strength F t 
that shifts the potential minimum from r — to the clus- 
ter surface r = Rt, i.e. [dVt(r) / dr] R = F t . From this the 
maximum charging of the cluster can be estimated to be 



Qt 



RtFt . 



(4) 



The depth of the minima in the energy of the electrons in 
Fig-fflhidicates directly the charging of the cluster. Hence, 
if Eq. (0J is true the depth and the early charging process 
in time should agree if the energy (j/-axis) in the panels is 
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scaled by F t . Indeed, for the three panels A-C belonging 
to T = 25 fs and T — 50 fs the temporal development of 
the electronic energy until the minimum is reached agrees 
very well after scaling (not shown). The considerations 
apply as long as the cluster has not expanded yet, i.e., 
Rt ~ Rq- For longer times the cluster expands and the 
scaling does not apply. 

3.5 Resonant energy absorption through collective 
electron motion 

If the pulse is long enough to allow for an expansion of 
the cluster during the illumination by laser light resonant 
absorption becomes possible. The mechanism behind res- 
onant absorption in a cluster of the present size is the 
match between the frequency of collective electron mo- 
tion within the cluster and the laser frequency. This leads 
to a characteristic phase lag of tt/2 between the oscilla- 
tion of the laser amplitude and the electronic response [2]. 
If one compares Fig. with Fig.Q]one sees that resonant 
absorption goes hand in gloves with a sharp decrease of 
the minimum electron energy (i.e. a fast charging of the 
cluster). This is clearly visible in all panels for T = 400 fs 
and in panel A for T = 200 fs. If the resonance occurs 
during the fall of the laser envelope (t > 0) it cannot be 
directly recognized in the electron energy (compare Fig.^ 
with Fig.0). 

3.6 Fine structure on the scale of the optical cycle 

Finally, we would like to draw attention to a peculiar be- 
havior in the periodic response of continuum electrons and 
bound electrons to the driving laser. Most intensity of 
such electrons can be found close to E = and close to 
the minimum of the electron energies, respectively. Firstly 
one sees that the spots of high electron density occur for 
bound and continuum electrons with a phase difference of 
half a period, cf. Fig. 03 Secondly, and may be even more 
puzzling on a first glance, the energy distribution of the 
electrons seems to "breeze" , with a narrow distribution at 
zero field and a wide distribution at maximum field for 
the bound electrons and exactly the opposite, although 
less pronounced, for the continuum electrons, cf. Fig. 01 

The behavior of the continuum electrons (Ei > 0) in 
Fig.|3| is easily understood. They have a certain spread 
of drift velocities p t or momenta as a result of the en- 
ergy absorption and the individual electron has an en- 
ergy of E t = (p, t + A(t)) 2 /2 with A(t) = zf drF(r) « 
zFt/uismut — zF t /u) cos(ujt — tt/2). In other words, since 
these electrons oscillate with a phase lag of it the veloc- 
ity or the momentum has a phase lag of tt/2. Hence, the 
spread of the p i is amplified due to the squaring in the 
energy at zero field, as can be seen from the peaks be- 
tween the white lines in the upper part of Fig. [3] and from 
the corresponding distribution of energies Ei > in Fig.|2] 
(black curve). 

The energies Ei of the bound electrons (Ei < in 
Fig.EJ show a similar spreading phenomenon, however 




Fig. 3. Magnification of the electron spectrum of series A at 
T = 50 fs from Fig.Q The white lines indicate times of maxi- 
mum field strength, the thick white line indicates the time of 
the cut through the spectrum shown in Fig.^] 

with maximum spread at maximum field. This can be 
understood if one takes into consideration that the many 
electrons in the cluster are an interacting system of charges 
which reacts very fast to forces, trying to reach an equi- 
librium position and thereby maintaining the spatial dis- 
tribution in the cluster. In a first approximation it can be 
assumed to be incompressible. This distribution is driven 
by the laser and has a vanishing phase lag since the eigen- 
frequency 0% = Q t / Rf of the cluster potential Eq. (j^J is 
larger than the laser frequency to. Hence, the total energy 
Ei of electron i is approximately 

*~ 2 2R t + 2 { r * + n\ ) + ^ ) \r j -r i \ 
«f + (5) 

To arrive at Eq. JSJ from Eq. @ we, firstly, replaced the 
ionic potential by the cluster potential Q and, secondly, 
neglected the small term quadratic in F. From Eq. (0) it 
becomes clear that by monitoring Ei we see essentially the 
coupling to the field, i.e., 

Ei = E® + r t F(t) (6) 

where Ej 1 ' is the energy which varies slowly in time. The 
spread in Ei results from the different positions fi of the 
electrons within the cluster, \¥i\ < Rt. Indeed, the spread 
is amplified through the laser field amplitude for maxi- 
mum field while the electron energy is focused to a small 
range at zero field, in accordance with the observation in 
Fig-El or more clearly in Fig.0] 



4 Summary 

We have analyzed the time-dependent distribution of elec- 
tron energies resulting from the illumination of a large 
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Fig. 4. Cut through the electron energy spectrum of series A 
at 50 fs in the optical cycle around maximum laser intensity 
I (time t = 0) indicated by a thick white line in Fig.|3] The 
black curve is for almost vanishing field strength, the red (grey) 
curve for the maximal field strength within the optical cycle. 
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cluster containing about 10 4 Xenon atoms with a short 
intense laser pulse at 780 nm wavelength. The spectrum 
shows a rich structure on the time scale of the optical 
laser period with many interesting details which would be 
worthwhile to probe in the future using attosecond tech- 
niques. This might be even possible in the near future since 
extended systems like a cluster show interesting electron 
dynamics on the 100 attosecond time scale while more 
tightly bound electrons in atoms or ions are too fast for 
direct probing with pulses of 10-100 attosecond length. 

Moreover, for now we do not have detailed knowledge 
about time resolved electron motion in clusters and such 
experiments would also reveal if our way of describing this 
dynamics theoretically is adequate. 
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